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WHAT IS CLAIMED IS : 

1. A method for the blind identification of sources within a system 
comprising P sources and N receivers, wherein the method comprises at 
least one step for the identification of the matrix of the direction vectors of 

5 the sources from the information proper to the direction vectors a p of the 
sources contained redundantly in the m=2q order circular statistics of the 
vector of the observations received by the N receivers. 

2. A method according to claim 1, wherein m = 2q order circular statistics 
10 are expressed according to a full-rank diagonal matrix of the 

autocumulants of the sources and a matrix representing the juxtaposition 
of the direction vectors of the sources as follows: 

C m ,x= A q Cm,sA q " (11) 

where £ m>s = diag([C^^ ,-.^££^£,1) is the full-rank diagonal matrix of the 
15 m = 2q order autocumulants C™' ,p s des sources, sized (P x F), and where 

Ag = [a® 9 ° ® a x a® {q ° ® a*] , sized (N<* x F) and assumed to be of 

full rank, represents the juxtaposition of the P column vectors 

[«/'-" ®o- 

3. A method according to one of the claims 1 and 2, comprising at least the 
20 following steps: 

0: the building, from the different observation vectors x(f), of an estimate 
C m r of the matrix of statistics Cm^ of the observations, 

1: the singular value decomposition of the matrix C mx , the deducing 
therefrom of an estimate P; A of the number of sources P and a square root 
25 c m., of C m.. t . for example in taking C™ = B s |L S | 1 / 2 where | . | designates 
the absolute value operator, where L s and E s are respectively the diagonal 
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matrix of the P 7 A greatest real eigenvalues (in terms of absolute value) of 
C m t and the matrix of the associated orthonormal eigenvectors; 

2: the extraction, from the matrix 6"* = [n T , r* T ] T , of the N matrix 
blocks T n : each block r n sized (N«i- D x F) being constituted by the Nta-v 
successive rows of starting from the "Nte- M (rt-l)+l"th row; 

3: the building of the N(N-l) matrices ® n \,n2 defined, for all 1 < m * n 2 < N, 
by 0ni,n2 = r n i # r n 2 where # designates the pseudo-inversion operator; 

4: the determining of the matrix Vsou resolving the problem of the joint 
diagonalization of the N(N-l) matrices 0 n i,n2; 

5: for each of the P columns b p of A\ A q) the extraction of the K = iVte-2) vectors 
b p (k) stacked beneath one another in the vector b p = [b P (l) T , b P (2) T , 

b p (K)T]T; 

6: the conversion of said column vectors b p (k) sized (N 2 x 1) into a matrix 
Bp(fc) sized [NxN] ; 

7: the joint singular value decomposition or joint diagonalization of the K = 
M?-2) matrices B p (k) in retaining therefrom, as an estimate of the column 
vectors of A, of the eigenvector common to the K matrices B p [k) associated 
with the highest eigenvalue (in terms of modulus); 

8: the repetition of the steps 5 to 7 for each of the P columns of A; A q for the 
estimation, without any particular order and plus or minus a phase, of the 
P direction vectors a p and therefore the estimation, plus or minus a unitary 
trivial matrix, of the mixture matrix A. 

4. A method according to one of the claims 1 to 3 wherein the number of 
sensors N is greater than or equal to the number of sources P and wherein 
the method comprises a step of extraction of the sources, consisting of the 
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application to the observations x(t) of a filter built by means of the estimate 
A\ A of A . 

5. A method according to one of the claims 2 to 4 wherein C m , x is equal to 
5 the matrix of quadricovariance Qx and wherein m = 4. 

6. A method according to one of the claims 2 to 4 wherein C m , x is equal to 
the matrix of hexacovariance Hx and wherein m = 6. 

10 7. A method according to one of the claims 1 to 6 comprising a step for the 
evaluation of the quality of the identification of the associated direction 
vector in using a criterion such that: 

D(A, A) = (ai, <x 2 , , a P ) 

where 

1 5 a p = min[d(a p , a,)] (17) 

and where d(u,v) is the pseudo-distance between the vectors u and v, such 
that : 

I H 2 

d( "' M = 1 - iJ^pv) (18 » 

20 8. A use of the method according to one of the claims 1 to 7 for a 
communications network. 
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9. A use of the method according to one of the claims 1 to 7 for goniometry 
using identified direction vectors. 



